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ABSTRACT 

The physical origin of high velocity cool gas seen in galactic winds remains unknown. 
Following I Wang| ( p~995a| ), we argue that radiative cooling in initially hot thermally-driven out¬ 
flows can produce fast neutral atomic and photoionized cool gas. The inevitability of adiabatic 
cooling from the flow’s initial 10^ — 10^ K temperature and the shape of the cooling function 
for T < 10^ K imply that outflows with hot gas mass-loss rate relative to star formation 

rate oi ^ ^ radiatively on scales ranging from the size of the energy 

injection region to tens of kpc. We highlight the (3 and star formation rate surface density 
dependence of the column density, emission measure, radiative efficiency, and velocity. At 
Tcooh the gas produces X-ray and then UV/optical line emission with a total power bounded 
by r\j 10 ^ if the flow is powered by steady-state star formation with luminosity L^. The 
wind is thermally unstable at Tcooh potentially leading to a multi-phase medium. Cooled winds 
decelerate signiflcantly in the extended gravitational potential of galaxies. The cool gas pre¬ 
cipitated from hot outflows may explain its prevalence in galactic halos. We forward a picture 
of winds whereby cool clouds are initially accelerated by the ram pressure of the hot flow, 
but are rapidly shredded by hydrodynamical instabilities, thereby increasing /3, seeding radia¬ 
tive and thermal instability, and cool gas rebirth. If the cooled wind shocks as it sweeps up 
the circumgalactic medium, its cooling time is short, thus depositing cool gas far out into the 
halo. Finally, conduction can dominate energy transport in low-/3 hot winds, leading to flatter 
temperature profiles than otherwise expected, potentially consistent with X-ray observations 
of some starbursts. 
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1 INTRODUCTION 


Galactic winds are crucial to the process of galaxy formation, eject¬ 


ing gas from galaxies ([Heckman et 
star formation, shaping the stellar 

al.[[1990[l, helping to regulate 
mass function and the mass- 

metallicity relation (Dekel & Silk 

T9861 IFinlator & Dave||2008| 

|Erb|2008|[Peeples & Shankar|2011 

1 , and enriching the intergalac- 

tic medium with metals ([Aguirre et al.[2001|[Oppenheimer & Dave[ 

[2006[ [2008[l. Hot winds heated by supernovae and stellar winds 


or active galactic nuclei may be present in both star-forming and 
passive galaxies at all redshifts. Quantifying their prevalence, dy¬ 
namical importance, and observational signatures is a key area in 
observational and theoretical studies of galactic winds. 

An important outstanding question is the nature of the cool 
clouds of molecular, neutral atomic, and ionized gas seen in blue- 
shifted absorption and emission in rapidly star-forming galaxies 
and starbursts at all redshifts. The prevailing picture is that these 
clouds are driven out of the host galaxy by ram pressure accelera¬ 
tion from a supernova-heated hot wind (e.g., |Veilleux et al.|2005| l, 
but it is unclear if clouds can be accelerated to the required veloc¬ 


ities before being shocked and shredded by hydrodynamical insta¬ 
bilities ( |Scannapieco & Bruggen|20 13) [Zhang et al.|2015| l. Addi¬ 
tional mechanisms for cold cloud acceleration have also been sug¬ 
gested, including momentum deposition by supemovae, the radia¬ 
tion pressure of starlight on dust grains ( Murray et al.|2005 |2010 


201 1| Hopkins et al.|2012|[^ang & Thompson|2012[|Krumholz &| 

Davis et al.|2014[ [Thompson & Krumholz|2014[ 


Thompson et al.|20T5)l, and cosmic rays (Breitschwerdt & Schmut- 


zler|1994||1999||Jubeigas et al.|2008|[Socrates et al.|200^ anasz| 


et al.|2013 


Here, we revisit the physics of the radiative cooling of hot, ini¬ 
tially adiabatic flows. We argue for a picture where the cool gas in 
galactic winds and halos precipitates directly from the hot wind 
as a result of radiative cooling, based on earlier work by |Wang| 
( |1995a|b| ), which was subsequently developed by |Efstathiou 
and|Silich et al. p003| 2004] l; |Tenorio-Tagle et al"l ~i 2003 


(200q 


2005 


|2007| t; IWiinsch et ~ ~(|201 Ijl, a nd suggested in ultra-luminous in¬ 
frared galaxies by Martin et^ ( |2015| ). 


In Section |2.1| we present an analytic discussion that shows 
that mass-loaded winds always cool on large scales. The key point 
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2 Thompson et al. 


is that adiabatic cooling lowers the gas temperature to a value 
where further decreasing the temperature leads to more rapid cool¬ 
ing (e.g., Fig. 0 In this way, the cooling timescale can become 
shorter than the dynamical time for expansion, the assumption of 
adiabaticity breaks down, and the flow rapidly cools, radiating its 
remaining thermal energy content. We provide simple scalings that 
give the cooling radius, radiative efficiency, emission measure, col¬ 
umn density, and velocity at the cooling radius. We further identify 
critical values of the star formation rate surface density and mass¬ 
loading rate for the cooling instability. Section [2^ shows when we 
expect conduction to dominate energy transport in hot winds, lead¬ 
ing to much flatter temperature profiles close to the central galaxy. 

In Section |2.3| we solve the general equations for a time- 
steady flow with an arbitrary cooling/heating function. The strongly 
cooled wind forms an extended region of cool photoionized out¬ 
flowing gas. We show that the extended gravitational potential of 
galaxies decelerates the cooled wind on ~ few — 100 kpc scales. 
Section [2^ discusses the linear instability of the cooling flow. We 
show that it is both convectively and thermally unstable, and that 
the latter should enhance density fluctuations by of order ~ 100 . 
In Sectionj^we discuss our results in terms of observations of sys¬ 
tems including NGC 253 and M82, ultra-luminous infrared galax¬ 
ies with high velocity line emission, high redshift rapidly star form¬ 
ing galaxies, and the absorption line systems both at z ~ 0 ob¬ 
served by COS-Halos ( |Werk et al.|2014| ), and at high-redshift ( |Stei-| 
|del et al.|2010| ). We note in Section |3.2| that the very strong depen¬ 
dence of the cooling radius and efficiency on the mass-loading rate 
[3 motivates a picture of galactic winds where the cool clouds are 
initially launched by ram pressure acceleration and are rapidly de¬ 
stroyed by hydrodynamical instabilities. This process seeds the hot 
flow with both a higher mass loading rate, that rapidly causes cool¬ 
ing on larger scales, and density fluctuations that grow rapidly as 
the thermal instability develops. In Section [33] we discuss the ori¬ 
gin of the cool gas in halos and we calculate the cooling time of 
the wind after it shocks on the surrounding circumgalactic medium 
while blowing a wind-driven bubble. Our conclusions are presented 
in Section |4] 


2 WIND COOLING 
2.1 Analytic Considerations 

One might expect that a hot spherical wind never cools radiatively 
since it rapidly expands to low density (n oc r“^). To see that it 
can, first note that adiabatic cooling of an ideal gas dictates that 
P/p^ = constant, which implies that T oc (7 = 5/3). 

As the super-heated flow initially expands from T 10 ® K, cool- 
ing is dominated by bremsstrahlung, the cooling time is tcooi oc 
T^/Vn oc the advection time is tadv ~ r/v oc r, so that 
the ratio tcooi/tadv oc grows slowly with radius. Thus, if the 
flow starts with tcooi/^adv > 1 , it remains adiabatic as long as 
bremsstrahlung cooling dominates cooling. 

The situation changes when the flow cools enough that metal 
lines dominate the gas emissivity at T < 10^ K. At these lower 
temperatures, which the flow inevitably reaches at sufficiently large 
distance as a result of adiabatic cooling, the shape of the cooling 
function changes, and can be approximated as (e.g., [Mac Low &| 
[McCray] 1988[|Draine|2011j ) 

A(T) = s/n^ - 3 X 10“^^ Tf (10^ < ^ < 10^ K)’ (1) 



Figure 1. Illustrative example of the cooling/heating function in photoion¬ 
ization equilibrium (PIE) with a metagalactic UV background ( jOppen-j 
jheimer & Schaye|2013||Wiersma et al.|2009) . The red line shows the contri¬ 
bution from H and He, whereas the purple line shows the contribution from 
metals. Bremsstrahlung dominates cooling at high temperatures. Heating 
dominates at low temperatures (the “spike” at ~ 10^ K shows where jej 
goes through zero). Equation ^ gives a simple approximation to the cool¬ 
ing function between 10^ < T < 10^ K for solar-metallicity gas in col- 
lisional ionization equilibrium (CIE) and over a range of gas densities in 
PIE. 


where T 7 = T/10^ K and where ^ is meant to capture the depen¬ 
dence of A on the metallicity of the gas: ^ = 1 for Solar compo¬ 
sition and ^ ~ 0.3 and ~ 3 for composition of 0.1 and 3 times 
Solar, respectively. In this regime tcooi oc and the ratio 

tcooi/tadv oc decreases rapidly with radius. Thus, even 

if the flow starts with tcooi/tadv > 1 , on large physical scales it 
has a chance to become radiative. Although the temperature scal¬ 
ing, normalization, and density dependence of the cooling function 
depends on the metallicity of the medium and whether or not it 
is collisionally ionized, photoionized, or in ionization equilibrium, 
the upturn in A(T) at low temperatures is unavoidable because of 
line cooling, even in a gas of primordial composition. An exam¬ 
ple of a representative cooling function for solar metallicity gas in 
photoionization equilibrium (PIE) with the meta-galactic UV back¬ 
ground is shown in Figure Note that the outflows powered by 
supernovae in rapidly star-forming galaxies may have super-solar 
metallicities, depending on the mass loading and metallicity of the 
surrounding ISM. More discussion of metallicity in the context of 
radiative cooling is provided in Section [O] 


As in 


Wang|(T995a|b| l (see also |Efstathiou| 2000 | jSilich et al.j 


|2003||2004| |, we adopt the wind model of Chevalier & Clegg < |1985| l 


(hereafter CC85) to evaluate the importance of radiative cooling. 
We deflne dimensionless parameters /3 and a such that 


Mhot = P Mi. 


( 2 ) 


and 


^hot = 3 X 10^^ ergs s ^ a M*, (3) 
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where the latter assumes that alO®^ ergs per core-collapse su- 
pernova (SN) is thermalized in the hot plasma, that there are 
10“^ SNe/yr per M©/yr of star formation, and that the star for¬ 
mation rate M* is measured in Mq/jt. For a spherical galaxy of 
size R, the CC85 model assumes a constant mass loading rate and 
energy deposition rate per unit volume for r < and that the 
flow is purely adiabatic and mass-conserving for r > R. Given the 
parameters a and /3, the solution then provides the radial run of 
pressure, density, and velocity as a function of radius at all r. The 
character of the solution is such that the Mach number Ad = 1 at 
R. Appendix [A|pro vides a discussion of the CC85 model. 


Using the CC85 wind model for a wind expanding into a flxed 
solid angle U, employing the cooling function given in equation 
Q- and assuming that r ^ R(v ^ constant), we And that 


, o in6 OL I ito.3 

tcooi ~ 3 X 10 yr-^ 




no 7 M*,io ^ 


and 


tadv ~ 1 X lO'^yr 



(5) 


function turns over at temperatures below ~ 10^ K (Fig. 0. For 
a ULIRG with M* ~ 100 Mq yr-\ ~ 0.3. For an LEG 
with larger R ~ kpc and M* ~ SOM© yr“^, /?min — 0.6. These 
scalings imply that radiative cooling of hot galactic winds may be 
prevalent both in local starbursts and in galaxies on the star forming 
main sequence at high redshift ( |Wuyts et al.|2011| ). 

Equation § can be re-written as a critical condition on (3, 
such that rcooi obtains at a radius r ^ R 

/ 1 1 \ 0.342 / ^ \ 0.2T0 

/3eHt ^ 0.73a°-^^° (R^T (-3^) ■ 

\ Tcool J 

( 8 ) 

The fact that /3crit is similar to /3min implies that the maximum 
value of the cooling radius is ~ 10 kpc for our nominal scaling 
parameters. Substituting, we And that the maximum value of the 
cooling radius is 

, • X 0.273 

max/x. n \ i r i 0.274 oO.728 f \ .r.. 

rcooi (/^ = /^min) ^ 15 kpc Q i?o .3 ( - 1 • (9) 

If the wind is mass-loaded sufficiently to cool radiatively (/3 > 
/3min), it does so on a scale smaller than . 


where i^o .3 = i^/0.3kpc, rio = r/lOkpc, M*,io = M^/IOM© 
yr“^, Q.A'k — U/ 47 r, and we have dropped the dependence on the 
mean molecular weight and Hydrogen fraction in tcooi (we assume 
— 1 throughout this paper). In equation (|^ for tcooi, note the im¬ 
portance of the star formation rate surface density E* = M^I'kR^, 
which varies among rap idly star-forming galaxies by more than 


four orders of magnitude { [Kennicutt 1998 1 . Systems with higher E 


have shorter flow cooling timescale ( Strickland & Heckman|2009T ). 

The choice of parameters for scaling the relations above (and 
below) is motivated in part b y M82, which has R 250 pc and 
M* 5 — 10M© yr“^ (e.g.. 


Rieke et al. 


1980 


Forster Schreiber 


|et al.|200T] |2003| l. In addition, from X-ray observations M82 has 
been infen-ed to have a ~ 1 and 13 ~ 0.2 — 0.5 ( [Strickland &| 
Heckman|2009 1. These values of Mi, and R are also roughly ap¬ 


propriate for the individual massive star-forming clumps seen in 
high redshift galaxies ( jGenzel et al.|201l| . We will see below that 
the propensity for radiative cooling and its radiative efficiency are 
strong functions of the host galaxy parameters. The applicability to 
individual systems and collections of systems is discussed in Sec¬ 
tion [3] 


Setting tcooi = tadv, we derive the critical radius beyond 
which the flow becomes radiative in ther ^ R limit: 

rcooi 4kpc ^ . ( 6 ) 

Note the very strong dependence on both a and [3. A small change 
in either can move rcooi signiflcantly. However, [3 cannot be ar¬ 
bitrarily small and still have rapid cooling. Because the cooling 
function has a peak and drops at temperatures below ~ 10^ K 
(e.g., Fig. 0 . and because tcooi/^adv decreases as r while 

the cooling function of equation Q is applicable, there is a mini¬ 
mum value of the mass loading parameter /3 such that rapid radia¬ 
tive cooling occurs anywhere in the wind: 

P 0.636 f Ro.3 U47r ^ 

Pmin ^ 0.64 a ^-— . (7) 

\Mi,,io e J 

For 13 < /3min, the flow remains adiabatic to arbitrarily large scales, 
and although tcooi /tadv decreases in the temperature range where 
equation Q is applicable, it never crosses unity before the cooling 


By setting rcooi = we derive the critical mass loading rate 
such that the cooling radius collapses to the scale of the launching 
radius R. For this calculation we use the general CC85 solution 
valid at r = it (eq. [^requires r ^ R). We And that 

4crit(rcooi = R)~ ■ dO) 


The expressions above can be inverted to yield a critical star for¬ 
mation rate surface density (E,^ = Mi^/nR^) for the flow to cool 
at any given radius r ^ R 


k, crit 


~ 3.5 


M® „0.267 7 10 kpc V a" ™ d 4 ,r 


yr kpc^ 


R 


V rcooi J e 


( 11 ) 


which is exceeded by starbursts in the local universe and by rapidly 
star-forming galaxies at high-redshift (e.g., |Wuyts et al.|2011j |. 

We can also estimate the critical star formation surface density 
required for the flow to cool at rcooi = R using equation ( p^ , 
which employs the full solution to the CC85 equations at r = R 
and not the limit r ^ R (see Appendix A; see also Appendix A of 
[Strickland & Heckman|2009| l. We And that 


E*, crit (rcooi = R) ^ 420 


M© 




( 12 ) 


yrkpc^ ^3.70 

is required for the cooling radius to reach rcooi — i^[^The lat- 
ter is achieved and exceeded in l ocal ULIRGs (e.g., [Condon et al. 


1991 [[Downes & Solomon|19^[Barcos-Munoz et al.|2015[[Sc 


ille et al.|2014[ |. As a result, radiatively cooling hot winds may be 
a natural explanation for the fast li ne emission seen in these and 


similar systems with very high E* (Soto et 


: al .|2012 


Martin et al. 


[2015| ). More generally, much of the star formation at high redshift 
occurs in systems that meet or exceed the value of E*, crit given in 
equation 0 . and we therefore expect the radiative cooling of their 
winds to be important. 


The fact that the minimum value of the mass loading required 


^ This value is somewhat higher than the estimate of E,^^ crit(^cooi = 
R) ~ 300 M© /yr / kpc^ obtained by setting rcooi = Ria equation 0> 
which is only strictly applicable in the rcooi ^ R limit. 
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4 Thompson et al. 


for radiative cooling to happen anywhere in the supersonic flow 
(/3min) and the maximum value such that rcooi ^ R (/^crit in 
eq.p^ differ only by a factor of a few suggests that the velocity of 
the gas when it cools is constrained. This follows from the fact that 
the asymptotic velocity of the hot flow is u ~ 10 ^ km/s. 

Substituting ^min into v we And a maximum velocity of the cooling 
flow of 

. • X 0.180 

fmax(rcooi) 1250 km/s ( ) , (13) 

\ MAtt iio.3 / 


which, for extreme ULIRG-lik e parameters (e.g., R — 5 0 pc and 


Barcos-Munoz et al. 


2015| ) reaches 


M* = 100 M© yr“ ; see, e.g. 

~ 2600 km/s, while for LBG-like parameters {R — Ikpc and 
M* = SOM© yr“^) is ~ 1200km/s. Similarly, substituting ^crit 
from equation into Voo we And a characteristic velocity for 
cooling on scale R\ 


^crit(^cooi — R^ — T20 km/s 


a ^ M*,io 

Q47r ^0.3 


(14) 


which reaches ~ 1250 km/s and ~ 700 km/s for ULIRG and LBG 
parameters, respectively. Although subject to the uncertainty and 
potential diversity in both a and Q in the above equations, these 
scalings imply that velocities inferred from absorption and emis¬ 
sion lines in galactic winds should have a well-deflned maximum. 
They also imply that if one selects galaxies of approximately the 
same size, then v should correlate with M* weakly. We return to 
the issue of velocity profiles in radiatively cooling winds in Section 
|2.3| where we show that the flow can decelerate significantly on 
~ 1 — 100 kpc scales because of gravity, and in Sectionj^ where we 
argue that because of potentially super-critical mass-loading (i.e., 
P > Pcvit) fountain flows with a broad range of velocities may 
develop. Even so, Umax (eq. is a well-deflned upper bound to 
the velocity of cooled gas in rapidly cooling hot winds that can be 
tested with observations, and Ucrit(/3 = Pcvit) (eqs. 


13 


and 


101 is a 


characteristic velocity for hot winds (or portions of hot winds) that 
are sufficiently mass loaded to cool at Ucooi = R and then escape 
to large scales. 


Both the scale of TcooI and the value of /3crit have additional 
direct consequences for observables. If the flow becomes radiative, 
then it can give up at most its total thermal energy content, which 
varies strongly as a function of radius. At r = i?, the Mach number 
M — 1 and the thermal and kinetic energy are comparable, with 
the total enthalpy flux exceeding the kinetic energy flux by a factor 
of 3. If /3crit(ucooi = R) obtains, we would therefore expect the 
total radiated luminosity as the matter cools to be approximately 

'^max ^ £hot/L* ~ 10“", (15) 

where 

L* - 10“ Lq M*,io (16) 

is the bolometric luminosity from star formation. This value of 77 
is consistent with observations of line emission seen in IMartin eti 
[n^ ( |2015| ) (e.g., their Figs. 14 & 15). However, for rcooi ^ R, the 
radiative efficiency of the cooling flow decreases rapidly because 
adiabatic cooling drives down the thermal content of the matter as 
the radius increases. In the limit rcooi ^ R, the Mach number can 
be approximated (for 7 = 5/3) as ^ 16^/^(r/i?)^/^, and the 


radiative efficiency of the flow can be approximated as 
1 


T] = 




5.0 X 10 


''"cool 

-3 P' 


a 


Lx 

~ 100^ 

Eh.ot 


.3.89 . 


,2.84 ( 

^ f^47r Ro.3 


( 17 ) 


where Lx is the luminosity radiated from the wind as it cools. The 
subscript X represents the fact that rcooi occurs at a temperature 
between 10®'^ and 10^ K and is thus in the X-ray or far-UV regime. 
In Section [TT] we present a more detailed discussion of the radia¬ 
tive efficiency in the context of the models computed in Section 
|2.3| Figure shows that the radiative power of the wind is not 
in fact dominated by the cooling region, but by the hotter, effec¬ 
tively adiabatic region at smaller radius (see also eq.|^. Even so, 
equation 0 gives a first estimate of how the radiative efficiency 
of the cooling material scales with the parameters of the problem. 
Note, though, that this estimate for the radiative efficiency is just 
the thermal energy of the flow at rcooi- Once the temperature drops 
to ~ 10^ — 10^ K, it will be subject to phoionization heating from 
the UV emission from the galaxy or the metagalactic UV back¬ 
ground. The material will then expand at approximately constant 
temperature, in photoionization equilibrium, and it will continue to 
radiate (this is the outer zone discussed in |Silich et al.|2003]|2004| l. 
We argue in Section [T^ that this cool outflowing material may form 
the cool gas seen in galactic halos at all redshifts. 


The strong dependencies of rcooi on the parameters of the sys¬ 
tem translate into even stronger dependencies for the gas density, 
column density (N = nr), and emission measure (EM = n^r) at 
the cooling radius: 


ricooi — 2.0 X 10 ^ cm 


^7.34 ^1-58 


/ M^,io 

\ f^47r 


2.58 


(18) 


iVcooi ^ 2.5 X 10 


.-2 




,10 


\LlA7rR0.3 


1.79 

, (19) 


EMcool 

-^ ~ 1.6 X 10 

pc cm“° 


^ 2.37 

/?5.37 

^0.3 


r M.,10 

\ LIa'k 


4.37 


( 20 ) 


The gas density at the cooling radius is important for assessing the 
applicability of CIE, the propensity of the medium to cool below its 
temperature in photoionization equilibrium, and in determining the 
final density of clouds precipitated out of the hot flow by thermal 
instability (see Section [2^. The column density is important for 
assessing the effective absorption optical depth of a given line for 
absorption-line studies of winds viewed either “down the barrel” 
or at large impact parameter with background quasars or galaxies. 
The column is also important for estimating the scattering optical 
depth of Lya radiation and other resonant transitions in escaping 
the expanding flow, and for studies of emission-line halos and their 
velocity profiles. Finally, the emission measure is useful for char¬ 
acterizing the expected surface brightness of the radiating gas in 
resolved systems and for emission line studies of warm/hot gas as 
it cools from ~ 10® '^ K to ~ 10^ - lO'^ K. 


Note that the scalings above imply strong diversity among the 
observed properties of systems; even a slight change in /3, a, LL, R, 
or M-), can lead to significant changes in ricooi, /Vcooi, and EMcooi- 
Conversely, both the range of /3 and the range of velocities for a 
cooling flow are tightly bounded. This narrow range of /3 will im¬ 
print itself on the column density and emission measures inferred 
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from observations. For example, substituting /3min into A^cooi, one 
finds a minimum column density at which cooling can occur of 

. • V 0.18 

AT Q /I inl8 -2 f \ 1 ^ 

iVcooi,min-3.4x 10 cm ^0.82- (21) 

If the limiting total gas column density one could detect in a given 
survey is 10 ^° cm“^, for example, this will significantly limit the 
mass loading of the systems that could in principle be seen and 
the velocities that material would be expected to reach. If galac¬ 
tic winds are better described by a wide range of (3 along different 
sightlines in a given system, as we argue may be the case in Section 
|3.1| then the fastest cool photoionized gas (eq. in a given ob¬ 
served system with have the minimum column density, and it will 
be given by equation ( [2T] ). Indeed, the strong diversity implied by 
our expressions for ncooi, ^cooi, and EMcooi may be partially mit¬ 
igated if there is a self-regulating mechanism that drives (3 /3crit, 
as one might expect if cool clouds from the host are incorporated 
into the hot fiow, as we discuss in Section 

Finally, note that although we focus here on radiative cooling 
on scales larger than the host galaxy size R, radiative cooling can 
also occur at r < where the mass and energy are injected (eqs.[^ 
and[^ Appendix [A| . |Tenorio-Tagle et al.| ( |2007] ) have studied this 
problem both analytically and numerically in the context of individ¬ 
ual super star clusters. Assuming uniform volumetric mass loading 
{fnhou g/s/cm^) and energy injection rates (ehot; ergs/s/cm^), and 
for a fixed value of ehot /they find a critical stagnation radius 
i^st that approaches r = 0 and R in the limit of low and high Ehot, 
respectively. Only material injected with Est < r < R participates 
in a supersonic outfiow. Material injected at 0 < r < Rst rapidly 
cools, does not enter the outflow, and perhaps gives rise to a new 
generation of star formation ( |Tenorio-Tagle et al.|2007||Palous et al.| 
|2014| ). As Rst approaches R, the total amount of material ejected 
in the hot outflow decreases, but the specific energy of the wind re¬ 
mains fixed because the ratio ehot /rrihot is fixed everywhere within 
the injection volume. 

In the models considered here, radiative cooling within the 
host galaxy will decrease the total amount of matter that partici¬ 
pates in the outflow, decreasing the gas density, and thus increasing 
the critical mass loading rate needed for cooling on large scales. 
However, the magnitude of this effect depends on the energy and 
mas injection profile. For example, in CC85 and |Tenorio-Tagle et| 

( |2007| ) for star clusters, ehot and rhhot are assumed to be con¬ 
stant throughout the volume, so that most of the matter and energy 
are injected very near r ^ R (Mhot ~ rfihotE^). Thus, Rst must 
be very near R to significantly change the total mass loading of the 
outflow. Adopting a perhaps more realistic model for galaxy scales 
so that TTihot oc ehot oc (Mhot(^) oc r; i.e., an isothermal 
sphere; see [Zhang et al.| 20 T^ , most of the matter and energy are 
still injected on scales of order ~ R. For more realistic disk-like 
geometries, or for extended mass and energy injection without a 
well-defined “edge” (e.g., an R beyond which ehot = rhhot = 0), 
more work is required to understand both the sonic point and the 
mass loading. 

Note that for the CC85 model and the cooling function 
adopted, the ratio tcooi/tadv(^ < R) is largest at r ~ E. Thus our 
estimate of /3crit {r — R) in equation |To| gives the limit such that 
cooling is important throughout the profile, from r = 0 outwards. 
This is then the limit that Rst R and we expect the model to 
break down completely because for {3 > ficvit the wind is radiative 
at all radii. 


2.2 Conduction 


As an aside, we note that conduction may strongly affect the tem¬ 
perature profiles of low-/3 outflows. In the energy equation, the net 
sources and sinks include both radiative losses and conduction. As 
we have shown above, radiative losses are important for /3 > 1 
(eqs.|7][^. The net local energy deposition rate from conduction 
is given by 

4ond = V-[/^(T)VT], (22) 

where h{T) ~ 5 x ergs cm“^ is the Spitzer con¬ 

duction coefficient. 


Setting aside the important uncertainty concerning the 
strength and orientation of the magnetic field and assuming the 
temperature profile follows T oc sls appropriate for the 

CC85 model on scales r > R, equation ^22) can be written as 
£cond = 44/^(T)T/9r^. The conduction timescale {nksT/Scond) 
is then 


4ond 1-3 X 10® yr 


M*,w /r \ io/3 


(23) 


The strong /3, a, and r dependences of this expression stem directly 
from the normalization of the temperature in an energy driven 
model (T (x a/P), the radial dependence of the adiabatic temper¬ 
ature gradient, and the strong scaling of the conductive flux with 
temperature (oc 


The conduction timescale tcond has a much steeper radial de¬ 
pendence than the advection timescale tadv — r jv. We thus expect 
conduction to dominate at small radii. Setting tcond = tadv (eq.[^, 
we derive a critical radius beyond which the flow is adiabatic, and 
below which the flow is dominated by conduction: 

reond ^ 160 pc RtC (24) 

The fact that rcond — R in this expression implies that /3 ~ 1 is 
the critical value of the mass loading parameter below which con¬ 
duction begins to dominate the temperature gradient at R. Setting 
i^cond = R, we find thaQ 

/3crit,cond ' . (25) 

V M^^io J 

That is, for /3 < /3crit, cond, Eond < tadv on scale R and we expect 
conduction to dominate advection out to a scale given by rcond - For 
P > Pcvit, cond, we expect conduction to have a minor effect on the 
resulting solutions. Comparing equation ( [^ with equation ([7]) we 
see that over virtually any range of P the standard CC85 model is 
invalid: for P > Pmin the flow is radiative, while for P < /3cond it 
is highly conductive, and Pmin ~ /^cond- 

A full exploration of the importance of conduction for \ow-P 
outflows is left for a future work, but here we note that a primary 
effect will be to flatten the temperature profile for r < rcond- In 
this regime we can ignore advection and radiative cooling, and thus 
the energy equation is integrated trivially to yield the power-law 


^ This estimate ignores the strong temperature gradient in the CC85 model 
at R, which is artificially steep — and would therefore further enhance 
the importance of conduction at R — because of the assumed sudden and 
artificial cutoff in the sources of energy and mass injection at R in the CC85 
model. 


© 0000 RAS, MNRAS 000, 000-000 




















6 Thompson et al. 


relation ( |Parker|1964] ) 

T(r)ocr“^/'^ < r < rcond), (26) 

which is significantly flatter than the expected for an adi¬ 

abatic flow. This may have significant consequences for observa¬ 
tions attempting to diagnose the dynamical significance of the hot 
flow with X-rays. For example, [Strickland et al.|(19971 ) find a flat 
temperature gradient potentially consistent with equation ( [26t in 
the M82 superwind on kpc scales. This may be due to conduction 
given the low value of /3 implied by the hard X-ray observations of 
[Strickland & Heckman| ( |2009) (they find /3 ~ 0.2 — 0.5). 

A caveat to our discussion of conduction is that the form of 
equation ( [^ is only valid if (1) the collisional mean free path of 
the electrons is smaller than the characteristic scale of the tempera¬ 
ture gradient Ae < dr/d\n T, (2) the electron-proton equilibration 
timescale is less than the advection time, and (3) the electron ther¬ 
mal speed (Ve) is larger than the flow speed v. Our estimates sug¬ 
gest that these are all satisfied for the /3 > 1 models we focus on in 
this work, but that (1) and (2) break down rapidly for /3 < 0.2 —0.5. 
Thus for the minimal mass loading expected from supernova driven 
winds composed only of supernova ejecta (/3 ~ 0.2) we expect the 
basic hydrodynamic assumption of CC85 to be invalid and a com¬ 
plete revaluation of the physics is in order. 


2.3 A More Complete Calculation 


To explore the dynamics of radiative cooling in high-(3 hot winds 
more fully, we solve the spherical steady-state wind equations with 
an arbitrary cooling/heating function and a simplified equation of 
state. Combining the Euler equations for mass, momentum, and 
energy conservation, we find that (e.g., [Earners & Cassinelli|19'^ 
see Appendix [B]) 


^ ^ ^ / Ve - 4cs \ D ( q \ 
dr 2r \ J CyT ) 


(27) 


dp 

dr 



\PD 

f « ) 

V ci - ) 

' r CvT ' 



(28) 


dr 


2 D 

rC^ 




+ 


Cvv 


2 2 
Crp — V 


(29) 


where D = {T/p) dPldT\^ = kETlm^, Cv = dEldT\^ = 
{3/2)kB/rup is the specific heat at constant volume, = 
dP/dp\^ and = dPjdplj. are the adiabatic and isothermal 
sound speeds, respectively, ^ = s/p is the net heating rate per gram 
(see Fig. [^, and Ve = 2GM{r)lr is the “local” escape velocity 
(see Appendix [B]). For the purposes of an exploration of parameter 
space, we use the cooling/heating tables assuming phoionization 
equilibrium from [Oppenheimer & Schaye[ ( |2013| t (see [Wiersma et[ 
|al.|2009| ), which employ the Haardt & Madau[ ( |2012| l meta-galactic 
photoionizing background. In PIE, ^ is a function of density and 
temperature. An example showing the contributions to cooling 
from H-i-He and metals separately is shown in Figure at a char¬ 
acteristic density. We neglect conductive energy transport (Section 
|2.2| ), which should only matter for the low-^^ ostensibly adiabatic 
models that are not the focus of the current study. Lastly, note that 
the wind equations above assume that the composition of the matter 
does not change as a function of radius, that continuity is obeyed, 
and that all emission is optically-thin. 


We focus on solutions that start adiabatic at r = with 
tcooi > tadv In this limit, the Mach number of the flow is A4(r = 


R) = 1, in accord with CC85 and the analysis presented in Wang 


( |1995a| ). We therefore solve the CC85 model for the flow at r = 
and then treat equations ([27])-([^ as an initial value problem and 
integrate to large radii using standard methods. For the purposes of 
illustrating rapid radiative cooling on large scales, we assume that 
the wind expands into vacuum. More discussion of a surrounding 
ambient medium is provided in Section 

Figure shows profiles of temperature, density, and velocity 
for fiducial models with M* = 10 M© yr“^, and R = 300 pc, 
close to parameters representative of starbursts like M82 and NGC 
253, or the massive star-forming clumps seen in high redshift 
rapidly star-forming galaxies. The red dotted lines show the so¬ 
lution without a gravitational potential (ve{r) = 0) and the blue 
solid lines show the solution assuming a singular isothermal sphere 
(Ve{r) = 2a) with velocity dispersion a = 200km s“^. All mod¬ 
els assume efficient thermalization a = 1, as is implied by obser¬ 
vations of M82 { [Strickland & Heckman[[2009[ ), and then a range 
of 13 from 0.2 to 1.6 in steps of 0.2. For large {3 we see the rapid 
decline in temperature from ~ 10^ K to less than 10^ K on scales 
of ~ 1 — 10 kpc. Here, photoionization heating balances cooling, 
as in the example cooling function shown in Figure The mate¬ 
rial then continues to flow outward at fairly high velocity. The right 
panels show that the rapidly cooling solutions have high column 
densities and a broad range of velocities, potentially explaining the 
broad emission line features seen in some systems (e.g., [Soto et al.| 
[2012|[Martin et al.|2015| Sectionj^. 

Note that the “adiabatic” models shown with 3 < 0.6 are 
meant for comparison with the more strongly radiative solutions. 
In fact, we expect these low-^^ solutions to be strongly affected by 


conduction (eq.|^ out to a radius given by equation (241 and will 
thus have much flatter inner temperature profiles (eq. 7R i. For the 
13 — 0.2 model, rcond — 1.8 kpc. A full solution with conduc¬ 
tion and radiative cooling requires a complete revision of the CC85 
model and is left for a future work. 


Returning to the importance of radiative cooling, we see from 
the lower left panel of Figure that the high-/3 models in an 
extended gravitational potential slow down significantly on large 
scales ( [Wang[[T995a| ). Formally, the /3 = 1 — 1.6 models spon¬ 
taneously develop a sonic point at ~ 80, 30, 10, and 7 kpc. This 
behavior follows from equation in the limit of strong cooling. 
Once the temperature drops to ~ 10^ K, the evolution becomes 
quasi-isothermal and we can approximate g ^ 0 at radii r > rcooi- 
Taking Cg ^ a and Cg v equation ^n\ shows that dvjdr < 0 
and the gas decelerates because of gravity. The critical temperature 
at which the flow starts to decelerate is just the virial temperature 
T ~ 3 X10® K cr|oo- Integrating from the cooling radius outwards, 
we derive the critical radius where the flow decelerates significantly 
after cooling: 


rsiow ^ rcooi exp 


~ 3.7 
.<^200 ^. 


(30) 


where rcooi is given by equation ([^, and where we have used the 
velocity at the cooling radius in order to obtain the numerator in 
the exponential. We see that small changes in /3, a, or a change 
rsiow dramatically. However, note that the detailed behavior of the 
gas during deceleration is sensitive to our assumption of a purely 
isothermal dark matter halo out to large scales. In numerical ex¬ 
periments we find that if we put in a more realistic NFW-like dark 
matter profile, the sonic point moves out to very large radius, or 
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Figure 2. Left Column: Temperature (top), density (middle), and velocity (bottom) as a function of radius for a models with M* = 10 M© yr“^, R = 0.3 kpc, 
and /3 = 0.2 — 1.6 in steps of 0.2, calculated without gravity (red dotted) and with an assumed isothermal gravitational potential with a = 200km s“^ 
(blue solid). Note the rapid decrease in temperature to less than 10^ K in the high /3 models at the cooling radius and the decrease in velocity on large scales 
in the extended gravi tational potential. All models assume the PIE with a meta-galactic UV background and solar metallicity jOppenheimer & Schaye|201^ 
[Wiersma et al.|200^ . Right Column: Temperature, velocity, and column density N = nr (cm“^) versus one another for the same models. The differential 
luminosity as a function of temperature for these models is shown in Figurej^ 
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does not occur within ~ 300 kpc for the parameters in Figure]^ In- 
creasing {3 slightly, however, can then cause strong deceleration to 
reappear. As is evident from equation ( [^ all changes in the param¬ 
eters of the wind and host galaxy will strongly affect this evolution. 

The models shown in Figure [^assume Solar metallicity Z = 
Zq and PIE in a |Haardt & M^au| ( |2012| ) metagalactic back¬ 
ground at redshift z — We experimented with changes to both. 
Lowering or raising the metallicity increases or decreases, respec¬ 
tively, the critical value of the minimum mass-loading required for 
strong cooling as one would expect (y^min; eq.0. For example, for 
Z = 0.2 Z© models with /3 = 1 and 2 are qualitatively similar 
to Z = Zo models with ^ ~ 0.6 and 1.4, respectively. Note that 
the shape of the cooling function allows for rapid radiative cooling 
even in primordial gas (H+He in Fig.0: a model with Z = 0 and 
/3 = 2 resembles the model with — 1 shown in Figure Con¬ 
versely, a model with Z = dZ© and — 1 has a cooling radius 
at the same location as the ^ — lA lines in Figure]^ Changes to 
the redshift for the PIE calculation have essentially no effect on the 
location of the cooling radius or the critical 13 for cooling, but they 
do change the post-cooling temperature in accordance with what 
one would expect for a harder metagalactic background. 

Over much of parameter space the cooling radius occurs at a 
sufficiently high density and temperature that employing a standard 
CIE cooling function instead of PIE does not change the cooling 
radius signihcantly. However, the post-cooling evolution is com¬ 
pletely different. In CIE, the gas expands adiabatically after cool¬ 
ing with T oc which differs qualitatively from the approx¬ 

imately isothermal evolution on large scales for the models shown 
in Figure 

Finally, note that the flow time to 100 kpc is > 100 Myr and 
we expect the assumption of steady-state conditions to be violated 
on this timescale, at least in starburst systems. We further expect 
the wind to expand into a local circumgalactic medium, sweeping 
up and shocking matter in a galaxy-sc ale wind-blown bubble (see 
Section |T^ . 


2.4 Thermal & Convective Instability 


[Balbus & Soker| ( p~989| ) (hereafter BS89) studied the case of ther¬ 
mal instability (TI) in a dynamically evolving background, which 
is the case applicable to TI in galactic winds. The key radial range 
where TI is important is where the cooling time is less than the flow 
time. Otherwise there is insufficient time for the instability to grow 
signihcantly. Under these conditions, both convection and isobaric 
TI can in principle grow. The ampliflcation of density fluctuations 
by convection is given by a factor ( [Balbus & Soker|19'8^ 



where 


N'^ir) 


2 rup / 1 dP \ ds 
5 ks \ p dr J dr ^ 


(32) 


N{r) is the Brunt-Vaisala frequency, which characterizes convec¬ 
tive instabilities, s is the entropy per unit mass, and the integral 
over time in equation U) is co-moving with the flow so that 
dt — dr /u[^ Evaluating this over the radial range where initially 


^ Anisotropic thermal conduction along magnetic fields lines modifies the 


hot galactic winds cool appreciably, we find that / |A^|dt is typi¬ 
cally <C 1. Thus there is insufficient time for convection to grow 
appreciably. This is fundamentally because the Brunt-Vaisala fre¬ 
quency introduced by radiative cooling in a wind is given roughly 
by Ntcool ^ where A4 is the Mach number of the flow and 
the flow only spends a time dt ~ tcooi at radii where cooling is 
significant. Note that at large radii once the flow comes into PIE 
and T ~ constant, dsjdr > 0 and the flow is convectively stable. 


If the cooling rate is assumed to scale as oc n^A(T), the am¬ 
plification of density fluctuations by isobaric thermal instability is 
given by a factor ( [Balbus & Soker|19'8^ 


Ati ~ exp 



■/ ainA^ 

, dt 

[ 5 ; 

V dlnTj 

tcooi _ 


(33) 


Evaluating this for the cases shown in Figure we find that 
Ati ~ 30 — 300 (for [3 = 1.0 — 1.6, respectively), with all of 
the contribution due to the small radial range where the flow cools 
abruptly from T ~ 10® K to ~ 10^ — 10^ K (see upper left panel 
of Fig. [^. Physically, as the flow cools from high to low temper¬ 
ature, each logarithmic interval in temperature (which corresponds 
to dt ~ tcooi) contributes an order unity contribution to the in¬ 
tegral in equation ( [^ , so that the net amplification is set by the 
few decades in temperature decrease during the rapid cooling of 
the wind. 


The fact that the amplification of density fluctuations by TI is 
limited to a factor of ~ 100 implies that for the medium to develop 
a multi-phase structure there must be initial fluctuations of finite 
amplitude present in the flow. Given the high initial sound speed in 
the thermally-driven wind, one might expect any such fluctuations 
to have been erased by the time rapid cooling sets in. However, the 
hot outflow quickly becomes supersonic and thus expands outwards 
more quickly than it can establish pressure equilibrium. In addition, 
the mixing of cool gas into the hot flow will likely imprint large 
inhomogeneities on the flow, which will seed both instabilities and 
rapid cooling (see Section [3^. A n analogous situation can be seen 
in Fig A3 of [Sharma et al.[ ( 2012[ ) where an inhomogeneous cooling 
inflow (rather than a wind) with tcooi < ^adv and initial density 
fluctuations generates significant multiphase structure. 


Another requirement for equation ( |^ to apply is that the fluc¬ 
tuations must be isobaric, rather than isochoric. This implies that 
there must be initial density fluctuations over a spatial scale 


L< 


Hcool 


M 


(34) 


where iTcooi is the width of the radial region where cooling is rapid 
(the temperature scale height) and M here is for the hot wind prior 
to the onset of very rapid cooling. Assuming this is satisfied, TI can 
in principle produce final densities of 




100 ricooi 



/ Tcooi 

VlO^K 


(35) 


where ricooi and Tcooi are the density and temperature of the hot 
flow at the cooling radius rcooi (see eq. [^. However, because 
all of the gas ultimately cools to PIE in our fiducial models at 
~ 10^ —10^ K there is nothing to maintain pressure equilibrium for 


convective instability criterion in dilute plasma such as the galactic winds 
of interest here (Quataert 2008). This does not change our conclusion that 
convection is relatively unimportant in galactic winds because the timescale 
argument applies even to conduction-mediated convection. 
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gas having n/ ^ ricooi- Gas at different density will have slightly 
different temperatures in PIE, which will allow some density fluc¬ 
tuations even if pressure equilibrium is satisfled. More importantly, 
some of the cold denser gas produced by TI may not have time to 
expand signiflcantly on the flow time. One way to see this is to note 
that if clouds occupy a fraction of 10“^ of the area of the sphere but 
have densities of ~ lOOncooi, they will carry the same mass flux as 
the original hot wind. With that Ailing fraction and a post cooling 
temperature of 10^ — 10^ K the gas will not have time to expand 
to All the sphere on the flow time. The proper physical picture may 
thus be one in which the cool clouds produced by TI in galactic 
winds are initially over-pressurized relative to their surroundings 
and then expand out at ~ 10s of km/s, eventually re-establishing 
pressure balance. Alternatively, magnetic and/or cosmic-ray pres¬ 
sure may be important in conflning the cool gas. 


3 DISCUSSION 

Two interconnected puzzles persist in the physics of cool gas in 
galactic winds and halos that may be solved by rapidly cooling, 
thermally-unstable hot winds of the type presented in Figure 
and Section The first is the presence of fast, but cool, outflows 
observed in emission and absorption in winds from rapidly star¬ 
forming galaxies across the universe. How is the molecular, neutral 
atomic, and ionized gas accelerated to the hundreds and even sev¬ 
eral thousands of km/s in some systems without being shocked, 
shredded, and incorporated into the hot flow? The second is the 
pervasive warm neutral atomic and ionized gas seen in the halos of 
galaxies at low and high redshift ( [Steidel et al.|2010] |Werk et ^ 
|2014| ). What is the origin of this gas, and how does it persist on 
large scales without rapidly accreting? 


3.1 Cool Gas in Winds 

The prevailing picture of the cool gas in winds is that it is entrained 
from the host ISM in the hot wind, and then accelerated from the 
galaxy by ram pressure (e.g., |Veilleux et al.|2005| l. In accord with 
this view, the soft extraplanar X-ray emission from wind galaxies is 
often interpreted as shocked interfaces between the cool clouds and 
the surrounding hot medium (e.g., | Strickland et al.|[2000l |2002| l. 
However, it is unclear how the cool gas clouds are accelerated to 
the velocities and large scales seen without being shredded by hy- 


drodynamical instabilities and incorporated into the hot flow ([Klein 

et al.|1994| [Cooper et al.|2009| Marcolini et al.|2005||McCourt et 

al. 

2015| [Scannapieco & Bruggen[|2015| [Schneider & Robertson 

20 ] 

14). 

Figure [^suggests an alternative interpretation of both the soft 


X-ray emission and the cool gas seen in galactic winds. In this pic¬ 
ture, a hot flow with /3 > 1 expands from the energy injection 
region and cools radiatively on large scales ( |Wang|1995a| ). We pro¬ 
pose that the high velocity cool clouds seen in absorption line trac¬ 
ers like NaD, Mg, and other optical and UV resonance lines are 
precipitated directly from the hot gas in a rapidly cooling flow. Fur¬ 
thermore, the observed extraplanar soft X-ray emission may also 
originate from cooling, and may not necessarily be dominated by 
interfaces and shocks as in the standard picture. 

Viewed “down the barrel” towards the central star-forming 


host, one would expect to see fast, predominantly blue-shifted ab¬ 
sorption lines of partially ionized gas and soft X-ray emission. 
Emission lines are also possible (see eq.|^. Indeed, recently [Mar- 1 
[tin et al.| ( |2015| l suggested that the very fast Lya emission seen from 
local ULIRGs could originate from direct radiative cooling of the 
hot gas in a CC85-like outflow ( jSoto & Martin||2012[ [Soto et ah] 
|2012| ). We And a critical surface density of star formation in the 
ULIRG range — S* ^ 500 M© yr“^ kpc“^ — for rapid cooling 
on a scale of a few hundred parsecs, given high thermalization ef- 
flciency (eq. [H- A particularly interesting prediction of the model 
is a direct connection between the velocity of the material and its 
mass loading. Higher /3 implies more rapid cooling closer to the 
source, but also smaller velocity since v oc in the CC85 

model (see eqs.[^and[^. For ULIRGs like those discussed in 
[Soto et aL] ( |2012| ) and [Martin et al.j ( |2015| ), these scalings imply 
rapid cooling and bright line emission at > 10^ km/s, as is ob¬ 
served. Although speculative, we note that for ULIRG-like param¬ 
eters and with f3 y /^crit (^q-fl^J? ^cooi and Ncooi (eqs.[T^and[T^ 
reach ~ 10^'^ cm“^ and respectively, implying that 

after thermal instability the material may self-shield and rapidly 
form molecules (n/ ~ lO^ncooi ~ 10^'^ cm“^; eq.|^, depend¬ 
ing on its dust content. The possibility that molecules may form 
directly in a rapidly cooling hot flow should be assessed by follow¬ 
ing grain sputtering, grain-gas cooling/heating, and line cooling be¬ 
low 10^ K. Much less speculative is the fast atomic line emission 
and absorption with maximum velocity given by equation |T3] ) that 
should be seen from rapidly star-forming systems if the emission 
measures and column densities are high enough, respectively. This 
picture of fast cool outflows may provide an explanation for the 


high velocities measured for objects like those from Tremonti et al. 
( [2007| ); [biamond-Stanic et ari ( [2012| ); [Sell et al.] ( |2014| ); [Geach et 


[^ ( [2014| t and should be tested with absorption line observations 
from a large variety of systems ( [Heckman et al.j 19^ . 


An observational prediction of the radiative wind picture ad¬ 
vocated here is that the observed X-rays arise from recombination, 
with a well-deflned temperature progression from hot to cool, pos¬ 
sibly with non-equilibrium ionization physics, and then a spatially 
extended region of photo-ionized gas together with the product of 
the non-linear evolution of the thermal instability in a supersonic 
background (Section \2A\ . Such a picture might help explain the 
strong spatial correspondence between the Ha and soft X-ray emis¬ 
sion in galactic winds ( [Strickland et al.|2004a|b| ). We highlight the 
fact that for a medium with constant mass flux cooling from high to 
low temperature, there is a particularly simple relationship between 
the total energy radiated per logarithmic interval in temperature and 
the mass flux, given by 


dL 

d\nT 


2 Trip 


(36) 


This expression follows from dropping the first term on the right 
hand side of equation ( [^ and takiim the limit that the medium 
is highly supersonic c^, c|-)[j Figure shows the differ¬ 

ential luminosity dLjdhiT along the temperature profile for the 
models in FigureThe expected power-law dLjdXv^T (x T in 
equation ( [^ is recovered for the high-^ models. We have numeri¬ 
cally verified that equation ( [^ is precisely satisfled in the cooling 
region during the precipitous drop in T. Thus, if the cooling region 


^ Note that in the subsonic limit /c^ = and the numerical 

factor on the right hand side of equation changes from 3/2 —)■ 5/2 
(e.g., |Fabian|1994| . 
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is identified in observations, equation provides a convenient 
expression for estimating the mass loss rate directly from X-ray 
observations. 


As a nearby observational example of a system that might 
showcase the onset of strong radiative cooling in a galactic wind, 
we consider the ~ 0.3 — 0.6 kpc long conical limb-brightened frus¬ 
tum defined in |Strickland et al.| ( |2000| ) along the minor axis to 
the south of NGC 253. Using the count rate in the energy range 
0.3 — IkeV (0.06 counts/s) or 0.3 — 2keV (0.07 counts/s), tem¬ 
perature (0.5 keV), and volume reported by | Strickland et al.| ( |2000| ) 
for this region (assuming a distance of D = 2.6 Mpc), one derives 
a total radiated fiux of ~ 6 — 10 x 10 “^^ ergs cm“^ s“^ and a hot 
gas density of order ~ 0.05 — 0.1 cm“^ for gas with unity filling 
factor in CIE with Solar abundances. Comparing to the gas pres¬ 
sure and density derived from Ha observations by [McCarthy et al.| 
( |1987| ) in the same region, [Strickland et al.j ( [200()| argue that the 
true hot gas density is ~ 0.1 — 0.3 cm“^. These densities imply a 
cooling timescale of ~ 3 — 9 Myr, whereas the advection time is of 
order Myr for gas at 500 km/s. In general, we find that our radiative 
wind solutions have tcooi/^adv ~3 — lOatr^l — 2R at tem¬ 
peratures of k^T ~ 0.5 keV and that the ratio decreases r apid ly at 
larger radius because A oc as described in Section [TJ] The 

NGC 253 observations therefore appear consistent with our mod¬ 
els in which radiative cooling becomes important at larger scales. 
We thus propose that the Ha and X-ray emission from this region 
and beyond may be the result of the non-linear development of the 
thermal instability described in Section [^4] The fact that the Ha- 
emitting clouds and the hot gas are in pressure equilibrium, and that 
the Ha flux increases beyond this region (see Fig. 5 of [ Strickland et[ 
[al.|2000| right panel) supports the notion that rcooi ~ 0.5 — 1 kpc. 
The extended multi-kpc soft X-ray and line-emitting halo would 
then be interpreted as the large-scale aftermath of strong radiative 
cooling and thermal instability on sub-kpc scales. 


Note that if this feature in NGC 253 were to represent a region 
of strong radiative cooling, the hot gas mass loss rate implied by 
equation ( [^ is quite small, of order ~ 0.01 — 0.03Mo/yr, which 
would imply a small value for [3 if M* ~ 5 M 0 /yr and a potentially 
very large wind velocity if a ~ 1. However, as estimated above, 
Aooi/^adv ~ 3 — 10 for an assumed velocity of 500 km/s, implying 
that equation does not apply. The estimated Mhot should be 
larger by the ratio tcooi//adv, which would increase the implied /3 
to ~ 0.1. The additional complications made clear by the beautiful 
X-ray image of the frustum from [Strickland et al.[ l |2004b[ ) is first 
the importance of the opening angle Q — a naive estimate indicates 
Q ~ 0.01 — 0.05 — and second the potential importance of non- 
spherical areal divergence as the fiow expands on ~ 0.5 — few kpc 
scales ( [Suchkov et al.[[T9^ [1996| l. The latter may be especially 
important as the hot wind expands laterally just above the disk and 
encounters the halo. 


Although the example of NGC 253 is illustrative, the details 
of what is observed will also depend on the photoionizing spec¬ 
trum seen by the gas, the mass-loading along any given sightline, 
and whether or not the flow maintains ionization equilibrium (e.g., 
[Breitschwerdt & Schmutzler|1994[p^^ . Clearly much more work 
is required to fully assess strong radiative cooling as an explanation 
for the Ha and soft X-ray halos of NGC 253, M82, and other local 
starbursts. In fact, [Strickland et ^ ( |2002| ) argued against the pos¬ 
sibility of radiative cooling for the origin of the cool gas in NGC 
253, but did so using the cooling times derived on multi-kpc scales, 
where we would argue r ^ rcooi- The story of M82 may be more 
complicated because the wind has likely overrun halo gas from 



Figure 3. dL/dlnT as a function of temperature for the wind models 
shown in Figure]^ The power-law form dLjdlnT oc T given by equa¬ 
tion in the limit of strong cooling is evident for the high-(3 solutions. 
Note that these models have Mi, = 10 Mq yr“^, which implies a bolo- 
metric luminosity of Lboi ~ 10^^'^ ergs s“^ and a corresponding energy 
injection rate of ^hot 10^^-^ ergs s“^, both of which are substantially 
off the top of the plot. Most of the injected energy goes to adiabatic losses 
rather than radiative cooling. 


the tidal interaction with M81, or because it simply does not fit 
with the model of a radiatively cooling wind discussed here. Work 
by [Hoopes et al. ([2003 jl indicates w eak cooling on large scales. 
In addition, [Strickland & Heckman[ ( [2009[ ) found that a ~ 1 and 
[3 ~ 0.2 — 0.5 in M82 based on the diffuse hard X-ray emission 
observed in its core. Although this value of {3 is nominally below 
our analytic requirement for strong radiative cooling on kpc scales 
for Solar metallicity gas, the wind may undergo distributed mass¬ 
loading and entrain material on the scale of R, as we discuss below 
in Section [3.2[ In this way, the small value of [3 measured in the 
core may be compatible with strong radiative cooling (higher (3) on 
somewhat larger scales. Much more work is required to assess the 
possibility of a radiative fiow in the many systems for which a care¬ 
ful comparison is possible (e.g., VV 114, Grimes et al.|2006[ Mrk 


231,[Veilleux et al.|2014|[Rupke & Veilleux|201 1[ and many more; 
[Heckman et al.[[1990[ l, but the example of NGC 253 implies that 
very high spatial resolution observations may be required to local¬ 
ize the cooling region in mass-loaded starburst winds with X-ray 
observations. 

An additional zeroth-order observational constraint on radia¬ 
tive hot winds comes from the total diffuse X-ray emission ob- 
serv ed fro m galaxies (as in [Zhang et al.[2014l l. As discussed in Sec¬ 
tion 


2.1 


the rapidly cooling wind can radiate at most F^hot power, 
which for steady-state star formation is ~ 1 % of the bolometric 


power of star formation (eq. 15 andll^, L* = eM*c , where e ~ 


7 X10“^ is an IMF-dependent constant. In practice, the radiative ef¬ 
ficiency is much lower than this 1%. For example, even the /3 = 1.4 
model shown in Figure|^radiates ~ 0.05^hot — 5 x 10“^ Li, dur¬ 
ing the cooling phase, whereas the /3 = 1 model radiates ~ 5 
times less. Figure summarizes these results. It also shows that 
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the dLjdhiT distribution is quite flat for low-^^ models and that 
most of the energy is radiated between 10® — 10^ ® K temperatures, 
where the flow is effectively adiabatic, even for high-/3. 

These numbers for the radiative power of cooling winds can 
then be compared with observations of the total (integrated) soft 
X-ray emission from rapidly star-forming galaxies. For example, 
[Mineo Q\~2^\2Qi\2) report the 0.5-2keV X-ray luminosity of star¬ 
forming galaxies to be Lx ^ 5 x 10 ^®ergs s“^, implying 
Lx!Li, ~ 5x10“®. Inan attempt to approximate the bolo- 
metric X-ray luminosities of their galaxies from 0.5-10keV, they 
And a number ~ 20 times higher. Lx jLi, ~ 10“^. These num¬ 
bers imply that rapid wind cooling could contribute signiflcantly to 
the total X-ray budget of rapidly star-forming galaxies in extreme 
cases, but even the most mass-loaded model in Figure does not 
violate an upper bound from observations. For flxed mass loading, 
system radius, and thermalization efficiency, we note that the ra¬ 
diative efficiency scales nearly linearly with the star formation rate 
(eq. [T^. A particularly interesting case for a more careful obser¬ 
vational comparison is NGC 6240 where the diffuse hard X-ray 
emission is larger ( |Wang et al.|2014| ). 


3.2 Distributed Mass Loading and Cloud Rebirth 

The fact that cool clouds are readily destroyed by KH instabilities 
motivates a scenario for distributed mass loading and cool cloud 
production in hot winds. In this picture, supernova ejecta first ther- 
malize on small scales and with a small hot gas mass loading fac¬ 
tor of /3 ~ 0.2 (as appropriate for the mass loss from supernovae 
alone) producing conditions for an initially metal-rich wind (Z ~ 
few times Solar for supernova ejecta, and a-element enriched). 
This picture is similar to that discussed in [Suchkov et al.|(l996| ). 
The initially low-jd outflow is very hot and fast since the velocity 
and temperature of the wind scale as u ~ 2000(a0.2/^)^/^ km/s 
and T ~ 9 X 10^(a0.2//3) K, respectively. This wind sweeps up, 
shocks, shreds, and incorporates cool gas on the scale of the host 
galaxy R, with four important consequences: (1) /3 increases, (2) T 
decreases, (3) v decreases, (4) Z decreases. Note that the first three 
of these effects decrease the cooling timescale relative to the advec- 
tion timescale on larger scales, increasing the chances for strong ra¬ 
diative cooling as in Figure]^ In particular, the increase in /3 drives 
the cooling radius rcooi rapidly inwards to small radii (eq. [^. In 
addition, the swept up clouds seed the flow with density perturba¬ 
tions, which may grow via the thermal instability at rcooi (eq.[^. 
In essence, the rapid destruction of cool clouds on small scales ini¬ 
tiates their subsequent formation on larger scales, an almost literal 
transmigration: self-induced cool cloud rebirth after death in an ear¬ 
lier form. 

This scenario for producing cool fast gas from hot winds is 
complicated by the very strong dependencies outlined in Sectionj^ 
and by the multi-dimensional and time-dependent character of real 
galactic winds. Indeed, rather than thinking of a galaxy as having a 
single a, /3, and Q, we would instead argue for a picture in which 
winds with a wide range of fd and a, subtending fractions of 47r, 
emerge from a given galaxy. Along some sightlines the amount of 
gas swept up is small, [3 stays close to ~ 0.2, a is high, and the very 
fast, hot, and adiabatic wind escapes to large scales. Along other 
sightlines the amount of material swept up is enough to quench the 
wind completely: the velocity decreases sufficiently and the flow 
becomes sufficiently radiative (/3 > ficviu eq.fTol, that the outflow 
turns into a mass-loaded fountain on kpc scales, t hus, in looking at 



Figure 4. Cooling timescale for the wind material as a function of radius 
assuming that it goes through a strong shock, as might be appropriate for the 
reverse shock in a galactic wind-blown bubble propagating in the circum- 
galactic medium for the models presented in Figure The sharp decrease 
at the end of the high-/3 profiles is from the shape of the profile just before 
the sonic point that develops in our solutions at large radii and reflects a 
breakdown of the time-steady assumption used throughout this work. 


the upper left panel of Figurej^one might imagine all (3^ emerging 
from a single system, and not just one; the observations will then 
be a convolution of these diverse outflows and their interactions. 

This picture has interesting consequences. It suggests that the 
most rapidly cooling material with the strongest mass-loading, and 
hence the highest column densities and emission measures (eqs.p^ 
[20| , that can in fact escape to large scales will have velocities 
correlated with the local escape velocity. In this way, the process 
of cool gas incorporation may be self-limiting: if too much gas is 
swept up from the host, the wind never emerges. This process may 
be related to the observation by |Martin| ( |2005| ) that the wind ve¬ 
locity and escape velocity are correlated. Also consistent with this 
picture, detailed absorption line profiles from wind galaxies in gen¬ 
eral imply a broad distribution of column densities and velocities, 
with the fastest material having the lowest column density and vice- 
vers a (e.g.,|Rupke et al.|2005a|b[ but there are notable exceptions, 
e.g., |Diamond-Stanic et al.||2012| l. A prediction of this picture is 
that the limiting [3 would be given in terms of observables by equa¬ 
tion with a weak dependence on host size and star formation 
rate at flxed a. 


3.3 The Fate of Cool Gas in Radiative Winds 


Studies of quasar absorption lines in the outskirts of foreground 
galaxies imply that massive halos of cool gas are prevalent at both 
high (z 1 — 2) and low (z ~ 0.2) redshift (Steidel et al.|2010[ 


ITumlinson et al.|2()TH|Rudie et al.|20T^|Werk et al.|2014| ). Mod- 

eling the absorbing medium as T ~ 10^ K gas in photoionization 
equilibrium with a background metagalactic radiation held implies 
a total mass in this cool component that is comparable to or larger 
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than the stellar mass of the central galaxy (see |Werk et al.|2014| for a 
detailed discussion). Covering fractions of low-ionization material 
are high out to ~ 100 — 200 kpc around both star-forming and pas¬ 
sive galaxies, while higher ionization OVI absorption is found more 
frequently around star-forming galaxies ( |Tumlinson et al.|[2011| |. 
The observed linewidths are hundreds (not thousands) of km/s. The 
densities inferred from photionization modeling, n ~ 10“^ cm“^, 
are far too low for the gas to be in pressure equilibrium with an am¬ 
bient medium at the halo virial temperature, which challenges an 
explanation based on thermal instability within a hydrostatic hot 
halo ^Mo & Miralda-Escudell 19961 [Mailer & Bullock||2004[ see 
Fig. 12 of jWerk et al.lM?] ). 

We suggest that the observed absorption may originate from 
rapidly cooling outflows of the type presented in this paper. The 
calculations in Figure show that mass-loaded flows driven from 
rapidly star-forming galaxies can radiatively cool on > few — 
20 kpc scale, that these flows can reach large distances (even as¬ 
suming an isothermal gravitational potential well out to 200 kpc), 
and that they slow down from their initial 500 — 1000 km/s ve¬ 
locities as a result of the extended gravitational potential. On these 
scales we expect a multi-phase medium as a result of the non-linear 
development of the thermal instability at rcooi (Section \1A) and 
the interaction with intervening gas (Section [3^ . It is notable that 
for 13 = Mhot /> 1 we expect more material to be ejected 
than to be formed into stars, and that it is in precisely this regime 
that cooling should be strong (eq. |7]). For the examples shown 
in Figure the column density of photoionized ~ 10^ K gas is 
~ 10^^ — 10^® cm“^ on ~ 100 kpc scales for the no-gravity (red 
dotted) calculations with /3 > 1. For metallicities of ~ 0.1 — IZe, 
corresponding metal-line column densities would be in the observ¬ 
able regime. The blue solid lines evolve to high column density 
as they “stop” at the sonic point, but here our time-steady calcu¬ 
lations break down and more detailed time-dependent calculations 
are necessary to explore the dynamics. 

The low velocities and relatively high columns on large scales 
appear consistent with the observations, but at the cost of assum¬ 
ing that all of the systems probed have recently had strong mass- 
loaded galactic winds. This is prima facie problematic for the pas¬ 
sive galaxy halos in jWerk et al.| ( |2014] ). However, a basic piece of 
wind physics might help explain why even early-type galaxies have 
wind-mass-loaded halos: galactic outflows do not expand into vac¬ 
uum, as assumed in the calculations of FigureIn fact, we expect 
a newly born, hot, fast galactic wind to expand into a pressurized 
circumgalactic medium (CGM), to sweep up and shock this mate¬ 
rial, forming a wind-driven bubble with a characteristic reverse and 
forward shock conflguration. A sketch is shown in Figure]^ In the 
simplest case, a contact discontinuity separates the shocked wind 
material from the swept up and shocked CGM. Simple arguments 
show that a wind-driven bubble would reach pressure equilibrium 
with a constant pressure CGM on ~ 40 — 80 kpc scales, assuming 
a CGM density of 10“^ cm“^ and a temperature equal to the virial 
temperature. In a hydrostatic gas halo with a falling CGM pressure 
profile, the bubble may well run away to large scales, depending on 
the CGM structure, the wind lifetime, and its energy and ram pres¬ 
sure. Importantly, wind-driven bubbles expand slowly compared to 
the free wind velocities that drive them. There is thus a chance that 
the remnant bubble formed by a cooling galactic wind would per¬ 
sist on multi-hundred kpc scales for Gyr timescales, thus allowing 
cool gas absorption around galaxies that are passive today. 

Such an explanation for cool halos at \ow-z in the sample of 


jWerk et ^ ( |2014| ) should be fully developed. An immediate ob¬ 
jection to such a model would be that as the cooled wind encoun¬ 
ters the reverse shock in the wind-driven bubble it will be shocked 
to high temperature, potentially shutting off cooling, and thus no 
low-temperature photoionized gas would be observed outside of 
the cooling wind region. Although our calculations do not address 
the issue of a wind-driven bubble directly, we can make a prelimi¬ 
nary estimate of relevance to the problem by asking what the tem¬ 
perature and density in the reverse shock region would be if the 
wind profiles from Figure [^shocked at any given radius. If we as¬ 
sume a strong shockwave at any position, we can then compute the 
post-shock cooling timescale. Figure shows the result. We take 
the post-shock density and temperature to be nRs(r) = 4nwind(^) 
and TRs(r) = (3/16)u^ind(^)^p/^s and we then calculate the 
cooling time 

fcooi, shock = A;RTRs/nRsA(TRs). (37) 

The models with and without gravity differ significantly. The high- 
/3 models without gravity (red dotted) have cooling times of order 
Gyr only on scales below ~ 10 kpc. As the material goes through 
the reverse shock, it is heated to high temperature where the cool¬ 
ing time is long. The high post-shock temperature follows directly 
from the high wind velocities that are maintained even after cooling 
in the no-gravity models. In addition, at large radius these models 
shock at very low wind density (see left middle panel Fig.[^, which 
further prevents cooling. 

In contrast, the models with gravity (blue solid) slow down 
significantly, increasing the density relative to the models with no 
gravity, and have lower post-shock temperatures because of the 
lower wind velocity. Both effects act to decrease the cooling time. 
We see that for /3 > 0.8, all of the models have post-shock cooling 
timescales less than the Hubble time out to scales of ~ 100 kpc. 
These conclusions would likely be strengthened by calculating the 
bubble evolution together with the wind evolution since the relative 
velocity of the reverse shock and the wind fluid will lead to lower 
post-shock temperatures and shorter cooling times. 

Hot galactic winds with /^ > 1 can thus populate the halos 
of galaxies with cool gas in two ways. First, the freely expanding 
wind cools on small scales, undergoes thermal instability, and then 
evolves on 100 Myr timescales to very large scales, as in Figurej^ 
Second, these cool decelerating winds will shock on the ambient 
CGM, driving a wind-driven bubble, and although the post-shock 
cooling times can be long, they are generically less than the Hubble 
time. We thus expect the post-shock region to cool radiatively. See 
Figure 

These pictures for the origin of cool gas in the halos of galax¬ 
ies do away with the notion of pressure equilibrium with a virial- 
ized hot halo. They circumvent the problem of maintaining a large 
mass of cool material well below the virial temperature in the CGM 
for many halo dynamical times. 


4 CONCLUSIONS 


Following Wang|(|1995a|bjl and a num ber of other works includ- 


ing |Efstathiou|2000||Silich et al.|2003 

|2004| [Tenorio-Tagle et al. 

20031120051120071 IWunsch et al.|2011 

[Breitschwerdt & Schmut- 


zler||1994| |1999| we investigate how initially hot adiabatic galac¬ 
tic winds cool radiatively on large scales. We present approximate 
scaling relations for the cooling radius, density, column density. 
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emission measure, and radiative efficiency (Section M- Each of 
these quantities depends strongly on the hot wind mass loading pa¬ 
rameter 13 = Mhot / . However, hot winds generically undergo 

radiative instability — in the sense that tcooi < tadv — for ^ > 0.5 
(eqs.[^and[^ with interesting and testable dependencies on the 
thermalization efficiency a, the scale of the star forming region R, 
and the star formation rate In particular, as RjMi, decreases, 
the 13 threshold for the cooling instability decreases. Among other 
implications, the relatively small range of [3 for strong radiative 
cooling implies a well defined maximum velocity given by equa¬ 
tion and a minimum column density given by equation 
with weak dependencies on the star formation rate and host galaxy 
radius. 

Using the time-steady wind equations for an arbitrary heating 
and cooling function, we present numerical solutions showing the 
temperature, density, velocity, and column density as a function of 
radius in Figure]^ both with and without an extended gravitational 
potential, and assuming PIE with a metagalactic UV background. 
Models with 13 > I (eq. undergo strong cooling on small scales 
with the temperature dropping by two orders of magnitude or more 
over a small fractional radial scale. The qualitative behavior of the 
solutions is largely independent of the metallicity of the gas and the 
character of the UV background. 

The cooled gas forms an extended photoionized region that 
can reach large scales, depending on the gravitational potential. 
Even in the absence of a surrounding medium, the radiatively 
cooled wind decelerates significantly, and in our models sponta¬ 
neously develops a sonic point, which signals the breakdown of 
our time-steady assumption. Equation ( [^ shows that the decel¬ 
eration profile is exponentially sensitive to the parameters of the 
wind and the velocity dispersion of the host galaxy. For the same 
values of /3 and the thermalization efficiency a, this strong depen¬ 
dence on a means that cooled winds will be completely unbound 
from dwarf galaxies, but strongly bound to massive cluster galax¬ 
ies. This physics may connect directly to the efficiency with which 
baryons are retained in halos of a given velocity dispersion, and 
thus inform discussions of the z = 0 stellar mass function and halo 
mass function. 

The primary predictions of this model for cool gas are the scal¬ 
ing relations derived in Section combined with the basic picture 
that the hot gas cools via recombination, with a well-defined tem¬ 
perature progression and differential luminosity given by equation 
( [^ and shown in Figure This picture of strong radiative cool¬ 
ing may be a natural explanation for the high velocity line emission 
seen in local ULIRGs ( |Soto et ar]|2012| ), as also emphasized by 
[Martin et al.| ( |2015] ). Equations (jll j) and show that the criti¬ 
cal star formation surface density required for the cooling radius to 
collapse to the scale of the star forming region is in the ULIRG 
range — S* ^ 10^ M© yr“^ — but that strong spatially ex¬ 
tended radiative cooling is also expected from LBG-like galaxies 
with E* > 10 M© yr“^. Strong radiative cooling may also provide 
an explanation for the exceptionally high velocity outfiows probed 


by Mg II absorption in blue post-starburst galaxies ( 

Tremonti etj 

|al.|2007[|Diamond-Stanic et al.|2012[|Sell et al.|2014 

1 . For exam- 


pie, the system J0905+5759 reaches a velocity of ~ 2500 km/s 
(Fig. 3 of jDiamond-Stanic et al.|20f^ , close to Umax in equation 
CD- Further tests of the model with “down the barrel” absorption 
line studies should be undertaken. The application of the model 
to local starburst galaxies like NGC 253 and M82 is tentative and 
should be more fully explored with detailed models compared to 



Figure 5. Sketch showing the basic evolutionary stages discussed and the 
three distinct cool gas components. The hot fast wind emerges from the host 
galaxy, accelerating cool clouds to small velocities before they are incorpo¬ 
rated into the flow. The hot wind with high mass-loading factor cools radia¬ 
tively at the cooling radius and undergoes thermal instability (Fig.[^. The 
cooled wind decelerates on large scales, forming an extended photoionized 
region. The wind then shocks on the inner edge of the wind-blown bubble 
driven into the surrounding CGM. The cooling timescale for the shocked 
wind can be signiflcantly shorter than the Hubble time (Fig.|^. The robust 
existence of multiple regions of cool gas in galactic winds and halos may 
explain its prevalence in observations of star-forming and passive galaxies. 
The sketch here resembles recent results from high-resolution galaxy for¬ 
mation studies (e.g., |Onorbe et al.|2015) . 

high-resolution X-ray observations on small scales, as we highlight 
in Section ITT] 

The rapid radiative cooling of nominally hot galactic winds 
may help explain observations at low and high redshift that show a 
significant amount of cool photoionized gas in the halos of galaxies 
(Section [T^ . A particularly appealing aspect of this picture is the 
mass budget itself. The total mass in the cool gas halos is inferred 
to be comparable to or larger than the stellar mass of the galactic 
host ( jWerk et ^|2014j ), which is in line with our prediction that 
j3 = Mhot/Mi, >0.5 — 1 (eq.|^ is required for strong radiative 
cooling. 

Importantly, there are two sources of cool gas in galaxy ha¬ 
los. The first is the cooled wind itself, the extended photoionized 
regions shown in Figure The second is in the reverse shock that 
must inevitably form as the wind sweeps up the surrounding cir- 
cumgalactic medium. Figure shows that the post-shock cooling 
timescales can be less than the Hubble time on 100 kpc scales if the 
wind decelerates. In this way the cool gas is reborn on larger scales 
and at later times as the bubble moves slowly outwards. 

Combining this picture of the halo with our discussion of dis¬ 
tributed mass-loading and ram pressure acceleration of cool clouds 
from Section |3.2| we arrive at a story of two cool cloud trans¬ 
migrations. Cool clouds within the host galaxy are first swept up 
by the ram pressure of the hot flow. These clouds are rapidly de¬ 
stroyed on small scales and at low velocities <10^ km/s by hy- 
drodynamical instabilities, as is well-documented in both idealized 
high-resolution numerical simulations, and in parameterized mod¬ 
els that take into account both the gravitational potential and the 
radial structure of the hot flow ( jScannapieco & Bruggen|[?015[ 
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[Zhang et al.|2015t . The hot flow, now with enhanced mass load- 
ing and density perturbations, cools radiatively on larger scales, 
forming an extended region of atomic and ionized gas moving at 
~ 10 ^ km/s (the first transmigration) and then decelerating in the 
host’s extended gravitational potential. The high velocity cool wind 
then shocks on the surrounding circumgalactic medium, destroy¬ 
ing the cool gas again, but because of the relatively short cooling 
times (Figure the cool gas is again reborn on Gyr timescales 
(the second transmigration). In this way we generically expect three 
spatially-separated components of cool gas in rapidly star-forming 
galaxies, as illustrated in Figure]^ 

Many additional issues are left to be investigated. A pressing 
theoretical issue with direct observational implications is the char¬ 
acter of the non-linear development of the thermal instability in a 
supersonic background. The medium is thermally unstable (Sec¬ 
tion [ 2 ^ and the initial density perturbations may grow by a factor 
of more than 100 at the cooling radius (eq. [^, but high resolu¬ 
tion multi-dimensional simulations are needed in order to resolve 
the question of whether or not the entire medium cools monolithi- 
cally or whether something akin to discrete cool clouds of a char¬ 
acteristic size precipitate out of the surrounding hot background. 
Our speculations on the fate of a cool wind-driven circumgalactic 
bubble also need to be confronted with both parameterized semi- 
analytic calculations and multi-dimensional numerical simulations 
with radiative cooling (e.g., |Onorbe et al.|2015[|Sarkar et al.|2015| l. 

A number of microphysical issues also need to be investigated. 
We have shown that conduction (Section [2^ can dominate the en¬ 
ergy transport in low-/3 outfiows (eq.[^, and our estimates suggest 
that the assumption of collisionality and thermal equilibrium for 
the ionized plasma break down. A consequence of this fact is that 
over essentially any range of /3 the standard CC85 model is invalid. 
Comparing equation ([25l with equation ([7|) we see that for ^ > 0.5 
the breakdown of CCH5 is a result of relative cooling, while for 
/3 < 0.5 the breakdown is a result of conduction. Analytic and 
numerical investigations are needed to understand the temperature 
profile in the low-/3 regime (eq.|^. These issues may well directly 
impact the interpretation of high-resolution X-ray observations of 
nearby starbursts like M82 ( [Strickland et al.|19^ . 
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APPENDIX A: THE CHEVALIER & CLEGG (1985) 
SOLUTION 


The CC85 model neglects gravity and radiative cooling, and as¬ 
sumes uniform volumetric energy and mass injection rates inside 
the spherical starburst volume r < R and adiabatic expansion for 
r ^ R. Assuming an ideal gas with adiabatic index 7, CC85 de¬ 
rived a self-similar solution for the Mach number AT (r) of the form 
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for r > Standard methods are used to solve these expressions 
for the Mach number at all radii. For 7 = 5 / 8 , equation gives 
the useful expression 

M{rf ~ (M^ > 1). (A3) 


The overall normalization for the density, pressure, and temperature 
are given by the total mass injection rate Mhot and the total energy 
injection rate ^hot (eqs. [^and|^. The velocity of the flow is given 
by 


2M^ 


F/hot 

M2 + 2/(7-1)J ^ 


(A4) 


the density follows from mass conservation (Mhot = Qr^pv), 
and the gas temperature can be computed from the adiabatic sound 
speed via M. For r ^ R, these are simply 
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to expand the radial pressure and specific energy derivatives, re¬ 
spectively, and then employing the identities 
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(BIO) 


and substituting into equations ( |B2| ) and ( |B3| ), we arrive at the three 
coupled ordinary differential equations given in equations ( [M] )- 
( [^ . Note that equation ( |B8| ) is identically zero for an ideal gas, 
but the general expressions are recorded here for completeness. 
Also note that these expressions neglect changes in the chemical 
composition or chemical potential of the matter. 


APPENDIX B: DERIVATION OF THE GENERAL WIND 
EQUATIONS 


The time-steady equations of hydrodynamics in spherical symme¬ 
try with optically-thin heating and cooling, and with only gravita¬ 
tional and pressure forces are 

M = Anr^pv = constant, (Bl) 


_ _ld^ _ GM{r) 
dr p dr ^2 ’ 


(B2) 


and 


de vP dp 
dr p‘^ dr 


(B3) 


where q is the net heating/cooling per gram, v is the velocity, p is 
the mass density, r is the radial coordinate, M(r) is the enclosed 
mass, e is the specific internal energy, and P is the pressure. Writing 
equation in differential form as 


d In p d In u 2 

dr dr r 

and employing the thermodynamic identities 
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and 
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^ See jZhang et al.| ( [20T^ for generalizations of equation to power-law 
energy- and mass-injection models for r < R. 
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